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Abstract. This paper develops and analyzes a semi-discrete and a fully dis- 
crete finite element method for a one-dimensional quasilinear parabolic sto- 
chastic partial differential equation (SPDE) which describes the stochastic 
mean curvature flow for planar curves of graphs. To circumvent the difficulty 
caused by the low spatial regularity of the SPDE solution, a regularization 
procedure is first proposed to approximate the SPDE, and an error estimate 
for the regularized problem is derived. A semi-discrete finite element method, 
and a space-time fully discrete method are then proposed to approximate the 
solution of the regularized SPDE problem. Strong convergence with rates are 
established for both, semi- and fully discrete methods. Computational ex- 
periments are provided to study the interplay of the geometric evolution and 
gradient type-noises. 



1. Introduction 

The mean curvature flow (MCF) refers to a one-parameter family of hypersur- 
faces {r t } t > C R d+1 which starts from a given initial surface r and evolves 
according to the geometric law 

V n (t,-) = H(t,-), 

where V n (t, ■) and H(t,-) denote respectively the normal velocity and the mean 
curvature of the hypersurface T t at time t. The MCF is the best known curvature- 
driven geometric flow which finds many applications in differential geometry, geo- 
metric measure theory, image processing and materials science and have been ex- 
tensively studied both analytically and numerically (cf. [TUl E3 EH] and the 
references therein). 

As a geometric problem, the MCF can be described using different formula- 
tions. Among them, we mention the classical parametric formulation |18| . Brakke's 
varifold formulation [2J, De Giorgi's barrier function formulation [SOU II] > the vari- 
ational formulation pQ, the level set formulation (5^1 HH [5], and the phase field 
formulation p~3l Qj)] . We remark that different formulations often lead to different 
solution concepts and also lead to developing different analytical (and numerical) 
concepts and techniques to analyze and approximate the MCF. However, all these 
formulations of the MCF give rise to difficult but interesting nonlinear geometric 
partial differential equations (PDEs), and the resolution of the MCF then depends 
on the solutions of these nonlinear geometric PDEs. One interesting feature of the 
MCF is the development of singularities, in particular singularities which may occur 
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in finite time, even when the initial hypersurface is smooth. The singularities may 
appear in different forms such as self-intersection, pinch-off, merging, and fatten- 
ing. To understand and characterize these singularities have been the focus of the 
analytical and numerical research on the MCF (cf. [51 HTJ 1 IT U [T5 ] fM\ 125 } 152] . and 
the references therein). 

For application problems, there is a great deal of interest to include stochastic 
effects, and to study the impact of special noises on regularities of solutions, as well 
as their long-time behaviors. The uncertainty may arise from various sources such 
as thermal fluctuation, impurities of the materials, and the intrinsic instabilities 
of the deterministic evolutions. In this paper we consider the following form of a 
stochastically perturbed mean curvature flow: 

(1-1) V n = H(t,-)+eW t , 

where W denotes a white in time noise, and e > is a constant. It is easy to 
check that (cf. [3TJ[T2]) the level set formulation of (1.1) is given by the following 
nonlinear parabolic stochastic partial differential equation (SPDE): 

(1-2) df = |V X ,/| div x / (ipT^) dt + e\V x ,f\ o dW t , 

where / = f(x',t) with x' — (x,Xd+i) denotes the level set function so that T t is 
represented by the zero level set of /, and 'o' refers to the Stratonovich interpreta- 
tion of the stochastic integral. Again, stochastic effects are modeled by a standard 
K- valued Wiener process W = {Wt; t > 0} which is defined on a given filtered 
probability space {J~u t > 0},P). 

In the case that / is a d-dimensional graph, that is, f(x',t) — x c i+\ — u(x,t), 
equation (1.2 1 reduces to 

y x u 



(1.3) du= y/l + |V^| 2 div z ( — V ;" \dt + ey/\ + \V x u\ 2 °dW t . 



To the best of our knowledge, a comprehensive PDE theory for the SPDE (1.3) is 
still missing in the literature. For the case d = 1, (1.3) reduces to the following 
one-dimensional nonlinear parabolic SPDE: 

(1.4) du= - a " — dt + e^l + \d x u\>odW t 

y/1 + \d x u\ 2 

= d x (axcta.n(d x u))dt + eyl + |9 x m| 2 o dWt- 

Here d x u stands for the derivative of u with respect to x. This Stratonovich SPDE 
can be equivalently converted into the following Ito SPDE: 

(1.5) du = \~dlu + (1 - f ) . % - 1 dt + ey/1 + \d x u\ 2 dW t 



d x u + (1 — ^-)&rcta,n(d x u)jdt + + \d x u\ 2 dW t . 



As is evident from ( 1.4 ), ( 1.5 1, the stochastic mean curvature flow 11.31) for d = 1 



may be interpreted as a gradient flow with multiplicative noise. Recently, Es-Sarhir 
and von Renesse |12j proved existence and uniqueness of (stochastically) strong 
solutions for ( |1.4[ ) by a variational method, based on the Lyapunov structure of the 
problem (cf. [121 property (H3)]) which replaces the standard coercivity assumption 
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(cf. [321 property (A)]). As is pointed out in [12], mild solutions for (1.4) may not 
be expected due to its quasilinear character. 

The primary goal of this paper is to develop and analyze by a variational method 
some semi-discrete and fully discrete finite element methods for approximating 
(with rates) the strong solution of the Ito form (1.5 1 of the stochastic MCF. The 
error analysis presented in this paper differs from most existing works on the nu- 
merical analysis of SPDEs, where mild solutions are mostly approximated with the 
help of corresponding discrete semi-groups (see [2D] and the references therein) . We 
also note that the error estimates derived in [T7] which hold for general quasilin- 
ear SPDEs do not apply to (1.5) because the structural assumptions, such as the 
coercivity assumption [17l cf. Assumption 2.1, (ii)] and the strong monotonicity as- 
sumption [171 cf- Assumption 2.2, (i)] fail to hold for (1.5|, and also the regularity 
assumptions [171 cf. Assumption 2.3] are not known to hold in the present case. 
In this paper, we use a variational approach similar to [171 [5J [7] to analyze the 
convergence of our finite clement methods. One main difficulty for approximating 
the strong solution of (1.5) with certain rates is caused by the low regularity of 
the solution. To circumvent this difficulty, we first regularize the SPDE 
adding an additional linear diffusion term Sd 2 u to the drift coefficient of (|1.5 
a consequence the related drift operator in (2.3 1 becomes strongly monotone, and 
the corresponding solution process u s is then i? 2 -valued in space. However, it is 
due to the 'gradient-type' noise that a relevant Holder estimate in the i/ 1 -norm for 
the solution u s seems not available, which is necessary to properly control time- 
discretization errors. In order to circumvent this problematic issue, we proceed first 
with the spatial discretization (3.1)- (3.2); we ma y then use an inverse finite element 
estimate, and the weaker Holder estimate (3.17) for the process u s h to control time- 
discretization errors. We remark that addressing space discretization errors first 
requires to efficiently cope with the limited regularity of Lagrange finite element 
functions in the context of required higher norm estimates, which is overcome by a 
perturbation argument (cf. Proposition 3.2). 

The remainder of this paper consists of three additional sections. In section[2]we 
first recall some relevant facts about the solution of ( |1.5[ ) from [12] ; we then present 
an analysis for the regularized problem. The main result of this section is to prove 
an error bound for u s — u in powers of S. In scction[3]we propose a semi-discrete (in 
space) and a fully discrete finite element method for the regularized equation (2.3) 
of the SPDE ( 1.5 1. The main result of this section is the strong L 2 -error estimate for 
the finite element solution. Finally, in section [4] we present several computational 
results to validate the theoretical error estimate, and to study relative effects due 
to geometric evolution and gradient-type noises. 



2. Preliminaries and error estimates for a PDE regularization 

The standard function and space notation will be adapted in this paper. For 
example, H 2 (I) denotes the Sobolev space W 2,2 (I) on the interval I = (0, 1), and 
H°(I) = L 2 (I). Let (v)j denote the L 2 -inner product on I. The triple (Q.,T,F) 
stands for a given probability space. For a random variable X, we denote by E[X] 
the expected value of X. 

We first quote the following existence and uniqueness result from [12] for the 
SPDE (1.5) with periodic boundary conditions. 
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Theorem 2.1. Suppose that u £ an d fix T > 0. Let e < y/2. There 

exists a unique strong solution to SPDE (1.4 1 with periodic boundary conditions 
and attaining the initial condition it(0) = uq, that is, there exists a unique H 1 - 
valued {J 7 t}te[0,T]-O'dapted process u = {u(t); t £ [0, T)} such that V-almost surely 



(2.1) (u(t),i P ) I = (uo,(p)j 



(d x ii, d x <p) j ds 
arctan(a x -u), d x (p) j 



ds 



Moreover, u satisfies for some C > independent of T > 0, 

(2.2) sup E[\\u(t)f m{1) ] <C. 

te[o,T] 

It is not clear if such a regularity can be improved from the analysis of [12] 
because of the difficulty caused by the gradient-type noise. In particular, ir- 
regularity in space, which would be desirable in order to derive some rates of con- 
vergence for finite element methods, seems not clear. To overcome this difficulty, 
we introduce the following simple rcgularization of (1.5 1: 

(2.3) +(!_£) * V 



dt + eJl + \d x u s \ 2 dWi 



To make this indirect approach successful, we need to address the well-posedness 
and regularity issue s for (2.31 and to estimate the difference between the strong 
solutions u s of ( 2.3 1 and u of (1.5|. 

Theorem 2.2. Suppose that Uq £ H 1 ^) and ||ito||#-im < Co, where Co > is 
independent of 8. Let e < \/2(l + 5). Then there exists a unique strong solution 
to SPDE (2.3) with periodic boundary conditions and initial condition u s (0) = Uq, 
that is, there exists a unique H 1 -valued {^Ft}te[o,T]-ddapted process u 5 = {u s (t); t £ 
[0, T]} such that there holds V-almost surely 



(2.4) (u\t),<p) I = (4, v >) I -{5+-) / (dvu'&tfjds 



ar 



ctan(d x u s ), d x ip) . ds 



l + \d x u s \ 2 ,(p ) dW s \/p£H 1 {L) Vte[0,T]. 



Moreover, u satisfies 

(2.5) sup E\\\\d x u s (t)\\l HI) 
te[o,T] lz 



+ SE 









\\dlu s {s)\\l, (I) ds 



< E 



^Molli 2 (/) 



2.1 



Proof. Existence of u s can be shown in the same way as done in Theorem 
(cf. [H]). To verify (2.5), we proceed formally and apply Ito's formula with /(•) = 
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|||^r • IIl 2 (j) t° ( a Galerkin approximation of) the solution u s to get 



where 



(*)IIW) 


+i 




\\\d x ul\ 


2 

mi) 


\\\d x u 5 o\ 






M t : 



V' 2 



2 ' -/n~'x u \\L 2 (I) 
J2 ft 



1- 



L 2 



c/.s 



2 Jo 



6U/1 + M 2 



L 2 



~ 2 d T U S -dlu S 2 



L 2 (7) 



ds + M t Vte[0,T], 



' (a^i + ia^ooi 2 ,^ 5 ^ 



is a martingale. Taking expectation yields 



E 



1, 



d x u 5 (t)\\l m + / [5\\dlu 



}2„,<5||2 
L 2 



dlu s 



y/l + \d x u*\< 



L 2J 



<E[^||^|| 2 2 ]. 



Hence, (2.5) hold. The proof is complete. 



□ 



Next, we shall derive an upper bound for the error u s — u as a low order power 
function of S. 

Theorem 2.3. Suppose that Uq = Uq. Let u and u s denote respectively the strong 
solutions of the initial-boundary value problems (1.5) and (2.3) as stated in Theo- 
rems 2.1 and 2.2. Then there holds the following error estimate: 

rT 



(2.6) sup E \\u d (t)-u(t)\\l HI) 
te[o,T] L 



+ 5E 









\\d x {u 5 (s)~u(s))\\l 2(I) ds 



< CTS. 



Proof. Let e 5 := u 5 - u. Subtracting ( pTT] ) from ( pQ| we get that P-a.s 
{e s {t)^) I = - [ \s(d x u,d x( p) I + (6+ e l)(d x e s ,d xV ) I 



ds + Mt 



+ (l — — ) (arctan(9 a; u' 5 ) — arctan^u), d x ip) 1 
for all ip £ H 1 ^) and t £ [0,T], with the martingale 

M t :=e^ + \d x u s \i - y/l + \d x vfi,v) I dW s . 
By Ito's formula we get 

(2.7) ||e tf (t)||i B(J) = -2jr t [*(S a! uAe a ) / +(*+^)||fl ie e a ||^ (f) 

e 2 

+ (l — — ) (arctan(9 a; u' 5 ) — arctan(3.j; it), d x e s ) j 



ds 



^\ + \8 x U^- y/l + \d x uf 



L 2 (7) 



ds 



2c 



I (yi + \d x u^ - vT+ |a^| 2 ,e 5 )^iy s 
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Taking expectations on both sides, and using the monotonicity property of the 
arctan function and the inequality (yl + : 



E[||e 5 (t)|| 2 L2(/) 



26 E 



< 6E 



\d x e 



<5||2 



Vl + y 2 ) <\x- y\ 2 yield 
ds < -26 E' 



(d x u,d x e ) I ds 



l^ U llL2(/) 



\d x e^ 



which and \2.2\ imply that 

6E 



E[\\e 5 {t)\\ 



mi)\ 



\d x e 



<5 11 2 



ds 



< 6E 
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ds 



\ d xu\\h(i) ds 



< (CT)6. 



The desired estimate (2.6| follows immediately. The proof is complete. 

3. Finite element methods 



□ 



In this section we propose a fully discrete finite element method to solve the reg- 
ularized SPDE (2.3) and to derive an error estimate for the finite element solution. 
This goal will be achieved in two steps. We first present and study a semi-discrete 
in space finite element method and then discretize it in time to obtain our fully 
discrete finite element method. 

3.1. Semi-discretization in space. Let = xq < x\ < ■ ■ ■ < xj+i — 1 be a 
quasiuniform partition of I = (0, 1). Define the finite element spaces 

V, h := {v h e C°(I); v h \ [x ., x . +l] e P r ([x j ,x j+1 ]), j = 0,l,--- ,j}, 

where P r ([xj, aJj+i] denotes the space of all polynomials of degree not exceeding 
r(> 0) on [xj,Xj + i]. Our semi-discrete finite element method for SPDE (2.3) is 
defined by seeking Uh(-,t,oS) : [0, T] x Q — > such that P-almost surely 

(3.1) (v£(t),v h ) I =(u e h (p) ) v h ) I -(5+£ 



(d x u s h , d x v h ) j ds 
(avctan(d x ul) , d x v h ) j ds 
+ €^^l + \d x u s h \ 2 ,v h )dW s Vv h €V r h Vt€[0,T], 



(3.2) 



u{(t,0)=u s h (t,l) Vie[0,T], 



where u s h (0) — P^itg, and P^ denotes the L 2 -projection operator from L 2 (I) to V^. 

To rewrite the above weak form in the equation form, we introduce the discrete 
(nonlinear) operator A s h : Vj? —> VJ? by 

e 2 

(3.3) (AfrWfaVhjj ■= (6+ —)(d x w h ,d x v h )j 

+ (l - — - ) (avctan(d x w h ),d x v h ) j VtOh, v h G V r . 



Then (3.1 1 can be equivalently written as 

(3.4) du{{t) = -A{u h (t) dt + eP h yi + \d xU {(t)\< 



dW, 
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Proposition 3.1. Fore < ^2(1 + 8), there is a unique solution u s h G C([0, T];L 2 (Q; Vj 1 )) 
to scheme (3.1). Moreover, there holds 



(3.5) sup E 

0<t<T 



i2 



SE 



\d x u s h (s)\r r T , ds 



i2 



< E 



4(o)" 2 



mi) 



e 2 T. 



Proof. Well-posedness of ( 3.4 1 follows from the standard theory for stochastic ODEs 
with Lipschitz drift and diffusion. To verify (3.5 1, applying Ito's formula to f(u s h ) = 



\ u h\\ 2 L 2 (i) an d using (3.4 1 we get 



(3.6) 



L\I) 



l4(0)H 2 r2 ,n-2 



L 2 (I 
t 



2c 



P^l + \d x u s h (s)\*f L 
Plsjl + \d x ufo)\\v$) ^ 



ds 



It follows from the definitions of A s h and P[ that 
(3.7) \\u{(t)\\h(i) < Il4(0)|l! 2(/ ) -(25 + e 2 ) j)\d x ui( s )\\ 2 L2{I} d S 

-(2-c 2 ) f (&rct&n(d x u s h (s)),d x u s h (s)) 



+ C< 



+ 2e 



t r 



ds 



1 + \\ d x u h(s)\\ L 2 {I) 



Then (3.5) follows from applying expectation to (3.7), and using the coercivity of 
arctan. The proof is complete. □ 

An a priori estimate for u s h in stronger norms is more difficult to obtain, which is 
due to low global smoothness and local nature of finite element functions. We shall 
derive some of these estimates in Proposition |3.2| using a perturbation argument 
after establishing error estimates for u s h . 

To derive error estimates for u s h , we introduce the elliptic if ^projection R r h : 
H l (I) ->■ Vf, i.e., for any w £ H\l), R r h w £ V r h is defined by 

(3.8) (d x [R r h w-w},d x v h ) I + (R r h w-w,v h ) I = Vv h £V r h . 
The following error bounds are well-known (cf. [5]), 

(3.9) ||«; - R r h w\\ L2(I) + h\\w - R r h w\\ H1{I) < Ch 2 \\w\\ H 2 {I) . 



Theorem 3.1. Let e < ^2(1 + 6). Then there holds 

r T 

(3.10) sup E \\u A {t)-u A h {t)\\ 2 ' 
te[o,T] 1 y ' 



+ SE 


11 







\d x [u {s)-u h {s)]\\ L2{I) 
< Ch 2 (l + 5- 2 ). 



ds 
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Proof. Let 



e d (t):=u\t)-u d h (t), r, d := u d (t) - R r h u d (t) , ? := R r h u d (t) - u 6 h (t) . 



Then e s = -q s + £ s . Subtracting (I3JJ) from fclM we obtain the following error 



equation which holds P-almost surely: 



.2 ft 



(3.11) (e 5 (t),v h ) I +(5+-) / (3^(8),^)^ 



^2 ft 



= - 1- 



arctan(9 a; u' 5 (s)) — arctan(d x u^(s)), c^Vh J 



f ( V /l + |9 x u 5 (s)| 2 -- V /l+ |^(s)| 2 ,^)^ s + (e 5 (0),« /l ) / 



for all w/j £ V 7 h . Substituting e s = rj S + and rearranging terms leads to 



(3.12) 



{d S (t),v h ) I + (5 + -) j (a^(a),^«fc)jdfl 

+ (l - y) y ^arctan(9 K u l5 (s)) - arctan(5 x u^(s)) , d x v h ^j ^ds 

dW s 



c 2 /■* 



/ (^(s),^) I rfs-(7 7 5 (t),v ft ) 7 +(e 5 (0),^) r 



Applying Ito's formula with f(£, s ) — M S \\ 2 L 2^y and using (3.121 and < 3.3 ) we 



obtain 



(3.13) \\t 5 (t)\\' L2{I) + (25 + e 2 ) J o \\d^ S ( S )f L2{I) ds 

+ (2 - e 2 ) ^ (arctan^-Rjyts)) - arctan(<3 x 4(s)), d x £ s (s)^ds 
= -(2 - e 2 ) ^ (arctan(d x i/(s)) - arctan (d x R r h u 5 (s)), a^s^ds 

+ e 2 /' Vi + |c^( s )| 2 - Vi + iaXOOl 2 " da 

Jo v v L U) 

+ 2e / + l^ 5 W| 2 - + I^^WP^W)^. 



(25 



<~) i ( V 5 (s),Z h ) I ds-2(r 1 5 (t),t 5 (t)) I +2(e s (0),Z 5 (t)) r 
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By the monotonicity of arctan, (3.9), (2.5), and the inequality (Vl + x 2 
\A + V 2 ) 2 < \ x — y\ 2 we have 



>o, 



' ^arctan(o x R^u 6 '(s)) — axctaxi(d x uf l (s)), d x £ s (s)J 
(2-e 2 )E J (arctan(d x u s (s)) - axctaxi(d x R r h u s (s)) , d x ^ 5 (s)^ ^ds 



< E 



\\d x e(s)\\h (I) +45-'\\d x u\s)-d x Rlu\s)\\l 2(I) 



Sf„\\\ 2 



>0 

5 - '* 



ds 



^4 E 



\\d^ 8 {s)\\ 2 LHI) ds +Ch z 5-\ 



E 



l + l^u^a)] 2 - Jl + |M(*)| 2 



L=>(J) 



ds 



< E 



"~ + 7 J / ll 5 *€ , (*)lli»(j) da 

«/ 



(»?'(*),£'(*)), 

(e'(0),^(t)) J 



)E 



< E 



+ C<T 1 E 
+ Ch 2 S- 2 , 



\\9xV 5 \\h(i) ds 



iU s (t)\\h {I) + \\v a (t)\\l HI) 



<1e 

- 4 



*AAI|2 

mi) 



< -E 
- 4 



■E 



llr(*)l 
i 



+ Ch 2 , 



\\e S (0)\\ 2 LHl) \ <- 4 U 5 (t)\\ 2 L2{I) +Ch 2 . 



Taking the expectation in (3.13) and using the above estimates then yields 



(3.14) sup E M 5 (t)f 



te[o,T] 



|^^(s)|| i2(J) ds 



< Ch 2 (l + 6- 



Finally, (3.10) follows from the triangle inequality, (3.9), and (3.14). The proof is 
complete. □ 



Remark 3.1. (a) Estimate (3.10) is optimal in the H 1 -norm, but suboptimal in 
the L 2 -norm. The suboptimal rate for the L 2 -error is caused by the stochastic effect, 



i.e., the second term on the right-hand side of (3.13), and it is also caused by the 
lack of the space-time regularity in L°°((0,T); H 2 (I)) for u s . 

(b) The proof still holds if the elliptic projection R r h is replaced by the L 2 - 
projection Pfc. 



We now use estimate (3.14) to derive some stronger norm estimates for u s h . To 
this end, we define the discrete Laplacian d 2 : Vj? — > Vj 1 by 

(3.15) (dlwh,v h )i = -(d x Wh,d x v h )i Vw h ,v h €Vf, 

and the L 2 -projection P£ : L 2 (I) -> V? by 



(PhW, v h )i = (w, v h )i Vv h €Vf- 
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Proposition 3.2. For e < i/2(l + 8) there hold the following estimates for the 
solution u s h of scheme (3.1 1: 



(3.16) sup E \\d x u{(t)\\l* ( i) 

0<t<T L 



-5E 



\\d 2 h u s h (s)\\l 2{I) ds <C{l + 5- 2 ), 



(3.17) E[\\u d h (t)-ui(s)\\ 2 L2{I) + -J \\d x [u d h (0 -ui(s)}\\i Hl) d( 

<c(i + s- 3 )\t- s \ yo<s<t<r. 

Proof. Notice that u s h = +R r h u s . By the ^stability of R r h , an inverse inequality, 
(f275j), and (pill) we get 



sup E 

te[o,r] 



\d*<(t)\\l H i) 



< 2 sup E 
te[o,T] 



d x R r h u s (t)\\h {I) 



< C sup E \\d x u 5 (t)\\ 2 LHl) 

t€[0,T] L 

< c(i + s- 2 ). 



C 

1? 



2 sup E 

te[o,T] 

sup E 

tG[0,T] 



d x C(t)\\l 2{I) 

2 



It follows from (3.15) and MM that 



II3/XHIW) = '{d x dlR r h w 1 d x Rlw) I 

= (d 2 h R r h w,d 2 x w) I + (w- R r h w,d 2 Rlw) I Vw e H 2 (I), 

and hence 

(3.18) \\d 2 R r h w\\ LHI) < \\dlw\lmi) + R r h w\\ LHI) < (1 + Ch 2 )\\w\\ Hm . 



By an inverse estimate, (3.18), and (3.14) we have 



E 



d 2 h u s h (s)\\h(i) ds 







< 2E 



(ch-^d x z\ S )\\i 2{I) +c\\dy(s)f L2{I) 

< C5-\l + 5- 2 ) + Ce\[ T \\d 2 x u s ( s )\\ 2 L2{I) ds 



< 2E 



ds 



which and (2.5| give the desired bound in (3.16). 

To show (3.17), we fix s > and apply Ito's formula to f(u s h ) — \\u 5 h {t) — 
u s h (s)\\ 2 L 2 ( j\ to get that for some {Ft] t e [s, T] }-martingale M t 



(3.19) \\u{{t)-u{{s)\\ 2 L2(I) 

= -(e 2 + 25) J (d x [u{{Q) ± u{{s)],d x [u{{Q - u{{s)]) ^ 

- (2 - e 2 ) f (arctan(^4(C)) ± arctan^u^)), d x [u{{£) - u{(s)}) d( 



\P^l + \d x [ul(()±ui(s)]\ 2 \\ 2 L2{I) d( + M t 
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By the L 2 -stability of P[, the triangle and Young's inequality, we can bound the 
last term above as follows: 



l 2 IL 2 (/) ^ 



; / \\P^l + \d x [ui(0±ui(s) 

J S 

'(\\d x [u s h (0 - u s h (s)]\\ L2(I) + ||1+ \d x ui(s)\\ 



< e 4 



< e 4 



L 2 (I) 



(i+s) \\d x [u d h (o-uUs)]\\i 2{I) dc 

+ e 2 (A + 6- 1 )(l + \\d x ui( S )\\l HI) ) \t- 



Also 



(£ 2 + 2<5) j (^<( S ),^(C)-t4(s)])/C 



< 



\\d*K(0 - <is)]\\i Hl) dC + (e 2 + 2SY6- 1 \t - s\\\d x u d h (s)\\t 2(I) , 



(2-e 2 ) J [axcUn(d x u s h (s)),d x {u 5 h {0-u d h (s)}) T dC 
6 '* 



< 



i 

2\2r-l 



\\d x [uUO - K(s)}\\ 2 LHI) dC + 4(2 - e^S- 1 \t - s\ 



Substituting the above estimates into (3.191 yields 

\\<(t)-u s h (s)f L2{I) + 6 - " 



\d x [ui(o-u s h (s)}fdc 



< 



C6- 1 ((2 - e 2 ) 2 + e 4 + 5 2 ) (l + \\d x u{{s)f L , {I) ) \t-s\ + M t 



Finally, (3.17) follows from applying the expectation to the above inequality and 
using (3.16|. □ 

3.2. Fully discrete finite element methods. Let t n = tit for n = 0, 1, • • • ,N 

be a uniform partition of [0,T] with r = T/N. Our fully discrete finite element 
method for SPDE (2.3) is defined by seeking an {Tt n \n = 0, l,--- , TVj-adapted 
V^-valued process {uJJ; n = 0, 1, • • ■ , N} such that such that P-almost surely 

(3.20) (t*J n+1 , v h ) I + r(S + t) (d x u S h n+ \d x v h ) I 

(5,n+lN 



+ r(l - — ) (arctan(9 s w^ n+ ),d x v h ) I 



(3.21) 



u h (°) = u h (!)> 



where AW n+1 := W(t n+1 ) - W{t n ) ~ 7V(0, r). 

We first establish the following stability estimate for w^'™- 

Proposition 3.3. Let e < ^2(1 + 5). For each n = 0, 1, ••• ,iV, i/iere is a 



valued discrete process {u S ^ n+1 ; < n < N — 1} which solves scheme (3.20) -(3.21 1. 
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Moreover, there holds 

S.n ||2 



(3.22) max E 

0<n<N 



h \\L 2 (I) 



25 



N 



n=0 



i2 

L»(J) 



< E 



<5,0||2 
X /i \\l 2 (I) 



e 2 T. 



Proof. The existence of solutions to scheme ( 3.20 1—( 3.21 ) for t, h > can be proved 
by Brouwer's fixed-point theorem, which uses the coercivity of the operator I + rA s h 
(see fl3~3"l». 

To show (3.22), we choose Vh = uf^ n+1 (uj) in (3.20) to find P-almost surely 



(3.23) 



5,n+l||2 



lL2(7) 
■2 



<5,n 1 1 2 
l h \\l 2 (I) 



1 1 1 5.n+l <5,n 1 1 2 

2lr/t IIl 2 (/) 



■(*+V)lM" +l " 2 



r 1 



<5,n+l <5,n 



-) (arctan^w,;™ ) } d x u^ 
AW n+1 . 



We compute 

(arctan(a 2; u^™ +1 ),a 2; u^" +1 ) / > 0, 

6,1 1 



^i + \d x u s h y,u s h n+1 



AW n+1 



^ 1 1 1 (5,n+l (5.n 1 1 2 

- 2ll u ' 1 ~ u h \\l 2 (i) 



1 + 1^%'' 



L 2 (I) 



| AW, 



n+l| 



By the tower property for expectations, there holds 



-E 



1 + l^xWft,'' 



L 2 (I) 



E[|AW n+1 | 2 |Ji„] 



= — tE 
2 



l^« 5 »"||! 2(/) 



such that we get 
1 



(3.24) 



E 



<5,ra+l II 2 



I £=> (J) 



-tE 



| a 5,n+l||2 



_i5,n 1 1 2 

IIl 2 (/) 

2 

L 2 (I) ~ 



tSE 



o <5.n 1 1 2 



la <5,n+l||2 



L 2 (I) 



L 2 (I) 



<e 2 r. 



After summation, we arrive at 



max E 1 1 1 u 

0<n<N L 



S,n 1 1 2 
h \\l 2 {I) 



N 

] +2<Jt^e[||cU£"| 



n=0 



2 

L2(7) 



< E 



5,0||2 
l h \\l 2 (I) 



So (3.22) holds. The proof is complete. 



e 2 T. 



□ 



Next, we derive an error bound for u s h (t n ) — u^ n . 
Theorem 3.2. There holds the following error estimate: 

N 



(3.25) sup E \\u a h (t n ) - u^ n \\ 2 2(n 

0<n<N L K ' 



+ 6E \Y1 T \\ d ^ U h( t n) ~ d x uS h l \\ 2 L 2 (I) 



< CT(1 



)h~ 2 T. 
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Proof. Let e s ' n := u{{t n ) - u^ 1 . It follows from §H§ that for all {t n ;n > 0} there 
holds P-almost surely 



(3.26) (u s h {t n+1 ),v h )j - (uft(tn),Uh)j 



^2 /■* 



(arctan^u^s)), d x v h ) jds 



(yjl + \d x u s h (s)\ 2 , » h ) dW. V« h G 



Subtracting (3.20) from (3.261 yields the following error equation: 



(3.27) (e s > n +\v h ) I -(a 5 > n ,v h ) I 

£ 2 rtn+1 

= -{S+—) J {d x u s h {s)~d x u } l n+1 ,d x v h ) I ds 

- (l - — ) / (arctan(d x u£(s)) - arctan(5 x u^" +1 ), 9 x u ft J ^ 



ds 



yi + |9^( S )P- y/l + dW.. 



Choosing = e 5, " +1 (cj) in (3.27) leads to P-almost surely 



(3.28) 2 [\\e s ^l 2{I) -\\e s ^l 2{I) }+^\ 

,2 



_<5,n+l <5,n||2 



I £2(7) 



+ (^ +5)r\\d x e^ +1 \\l H i ) 



(~ + S) f " +1 (d x ui(s) - d x u s h (t n+1 ),d x e 5 ' n+1 ) 



ds 



2 ft n+1 



-(!-—) y ^axctan(9 x itft(s)) ± arctan(<9 x u h (t n+ i)) 



arctan(^4 n+1 ),9 :E e^™ +1 )^ 



We now bound each term on the right-hand side. First, since E[AW n +i] = 0, 
by Ito's isometry, the inequality (-\A + x 2 — \J\ + y 2 ) < (x — y) 2 , and the inverse 
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f*" +1 (yJl + \d x u s h ( S )\* - y/l + fatfiy, e S - 
< E [i|| e ^" +1 - e 5 >l| 2(J) ] + j J' n+1 \\d x [u s h (s) - ± u S h (t, 



inequality we get 
E 



: PROHL 



,n+l ±e 5,„\ ^ 



Hi 2 (/) ^ s 



< -E 
~ 2 



~ 2 



jl^" +1 -e 4 * B ||i. w ]+E[(^ + |)r||^»||i 

IsXOO-u^tn)]!!! 



f 2 9 / ' 1 



2 T i 

|j e ^+i_ e ^||| 2(/) 



L 2 (/) ^ s 



(3.29) 

+ C(l + 5- 1 )h- 2 E 
An elementary calculation 



,e 2 5, 
+ (2 + 2-) TE 



ll^e^"||l 3(J) 



^" +I ||4( S )-4(t n )||2 2(J) d S . 
and an application of an inverse inequality yield 



j + 5) J^ +1 (d x u s h (s) - d x u s h (t n+1 ), d x e s > n+1 



ds 



< S -r ||a x ^» +1 || 2 L2(/) + 2( i± 5) ! \\d x [u{(s) - u{{t n+1 )]\\h {1) ds 

(3.30) < 5 -r Re^ +1 || 2 2(/) + (e 2 + 26) 2 S- 1 h^ \\u 5 h (s) - u s h (t n )f L2{I) ds. 



By the monotonicity of arctan we get 

,2 /•*„+! . 



£ 2 /•*„+! 

- (l - — ) y (arctan^u^s)) ± axcttm(d x u s h (t n+ i)) 

- arctan(9 2: u^ n+1 ), 9 a e' 5 '" +1 ^ 



-2 /■« 



<|r||S B e^+ 1 ||| 



- arctan(aX' n+1 ), d x e S ' n+1 j ds 

(axctan(d x u s h (s)) - axctim(d x u 5 h (t n+ i)) , 
2. e 2 z-*^ 



ll^[«h(«)-«A(*n+l)]lli»(/) ^ 

(3.31) < ||9 x e^ +1 ||| 3(J) + (2 - e 2 ) 2 5^h' 2 \\u s h (s) - <4(*»)ll W) 



Finally, substituting the above estimates into ( 3.28 ) , 
1, and using (3.17) and the fact that e S '° = we get 

jE[r^||9 x e^||! 



summing over n = 0, 1, 2, • • ■ 

e s '° = we get 



sup E 

0<n<N 



)\e S ' n \\l Hl) 



n=0 



<C(l + ^ 1 )/ l - 2 ^ / E sup ||4(s)-4(*n)ll! 2 (/) 



ds 



< CT(1 + S- 1 ) 2 h~ 2 r 
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which infers (3.25). The proof is complete. 



□ 



We conclude this section by stating the following error estimates for the fully 
discrete finite element solution u^ n as an approximation to the solution of the 
original mean curvature flow equation (1.5) 

Theorem 3.3. Let u and ui 



denote respectively the solutions of SPDE (1.5) and 
scheme (3.20) -(3.21 1. Under assumptions of Theorems 2.1, 3.1, and 3.2, there 
holds the following error estimate: 



(3.32) sup E — 

0<n<iV L 

< CTS + C(l + 



L 2 (I) 



N 

d x u(t n ) - d x u,; n \\ L2iI) 



n=0 



CT(1 



5- 2 )h- 2 r. 



Inequality (3.32) follows immediately from Theorems 2.3 3.1 and 3.2 and an 
application of the triangle inequality. 

Remark 3.2. We note that the main reason to have a restrictive coupling between 
numerical parameters in (3.32) is due to the lack of Holder continuity (in time) 
estimate for d x u s h in I? -norm. On the other hand, it can be shown that, under a 
stronger regularity assumption, the estimate (3.32) can be improved to 



(3.33) sup E \\u(t r , 

0<n<N L 



S,n 1 1 2 
U h IIl 2 (/) 

N 

5E["y^ T\\d x u(t n ) - d x i4' n ||^ 2(J) 

n=0 



< C(h 2 + r) 



This is because we no longer need to use the inverse inequality to get (3.29) -(3.31 1, 
and (3.33) can be obtained by starting with a control of the time discretization first. 



4. Numerical experiments 

In this section we shall first present some numerical experiments to gauge the 
performance of the proposed fully discrete finite element method and to examine 
the effect of the noise for long-time dynamics of the stochastic MCF of planar 
graphs, and we then present a numerical study of the stochastic MCF driven by 
both colored and space-time white noises where no theoretical result is known so 
far in the literature. 

4.1. Verifying the rate of convergence of time discretization. To verify the 
rate of convergence of the time discretization obtained in Theorem |3.3| in this 
first test we use the following parameters e — 1, 6 — 10~ 5 , and T = 0.1. In 
order to computationally generate a driving reference R-valued Wiener process, we 
use the smaller time step r = 10~ 5 . The initial condition is set to be u (x) — 
sin(7ra;). To calculate the rate, we compute the solution u.'" for varying r = 
0.0005,0.001,0.002,0.004. We take 500 stochastic samples at each time step t n 
in order to compute the expected values of the I/°°(L 2 )-norni of the error. The 
computed errors along with the computed convergence rates are exhibited in Table 
[T] The numerical results confirm the theoretical result of Theorem 13.21 
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Expected values of error 


order of convergence 


dt=0.004 


0.41965657 




dt=0.002 


0.27206448 


0.62526 


dt=0.001 


0.18136210 


0.58508 


dt=0.0005 


0.12373884 


0.55157 



Table 1. Computed time discretization errors and convergence rates. 



4.2. Dynamics of the stochastic MCF. We shall perform several numerical 
tests to demonstrate the dynamics of the stochastic MCF with different magnitudes 
of noise (i.e., different sizes of the parameter e). 

Figure [l] shows the surface plots of the computed solution u S ^ n at one stochastic 
sample over the space-time domains (0, 1) x (0, 0.1) (left) and (0, 1) x (0, 2 8 x 10~ 5 ) 
(right) with the initial value u (x) = sin(-7ra;) and the noise intensity parameter 
e = 0.1. The test shows that the solution converges to a steady state solution at 
the end. 




Figure 1 . Surface plots of computed solution at a fixed stochastic 
sample on the space time domains (0, 1) x (0, 0.1) (left) and (0, 1) x 
(0,2 s x 10~ 5 ) (right). uq(x) = sin(7ra) and e = 0.1. 



Figures [2][4] are the counterparts of Figure [T] with noise intensit y pa rameter 



e = 1, v2, 5, respectively. We note that the error estimate of Theorem 3.3 does not 
apply to the latter case because the condition e < y/2(l + S) is violated. However, 
the computation result suggests that the stochastic MCF also converges to the 
steady state solution at the end although the paths to reach the steady state are 
different for different noise intensity parameter e. 

We then repeat the above four tests after replacing the smooth initial function 
uq by the following non-smooth initial function: 



(4.1) u Q (x) = 



10a:, if x < 0.25, 

_ J 5 -10a;, if 0.25 < x < 0.5, 

10a: - 5, if 0.5 < x < 0.75, 

10 -10a;, if 0.75 < a: < 1. 
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Figure 2. Surface plots of computed solution at a fixed stochastic 
sample on the space time domains (0, 1) x (0, 0.1) (left) and (0, 1) x 
(0, 2 s x 10" 5 ) (right). u Q (x) = sm(nx) and e = 1. 




Figure 3. Surface plots of computed solution at a fixed stochastic 
sample on the space time domains (0, 1) x (0, 0.1) (left) and (0, 1) x 
(0,2 s x 10~ 5 ) (right). u (x) = sm{irx) and e = \?2. 




Figure 4. Surface plots of computed solution at a fixed stochastic 
sample on the space time domains (0, 1) x (0, 0.1) (left) and (0, 1) x 
(0, 2 s x 10~ 5 ) (right). u Q (x) = sm(irx) and e = 5. 

The surface plots of the computed solutions are shown in Figures [5}|8j respectively. 
Again, the numerical results suggest that the solution of the stochastic MCF con- 
verges to the steady state solution at the end although the paths to reach the 
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steady state are different for different noise intensity parameter e. As expected, the 
geometric evolution dominates for small e, but the noise dominates the geometric 
evolution for large e. 





Figure 5. Surface plots of computed solution at a fixed stochastic 
sample on the space time domains (0, 1) x (0, 0.1) (left) and (0, 1) x 
(0, 2 8 x 10" 5 ) (right). u Q is given in (gjl and e = 0.1. 





Figure 6. Surface plots of computed solution at a fixed stochastic 
sample on the space time domains (0, 1) x (0, 0.1) (left) and (0, 1) x 
(0,2 s x 10~ 5 ) (right). u is given in (lO) and e = 1. 



4.3. Verifying energy dissipation. It follows from (2.5) that the "energy" J(t) : 



|E[||9 a; M' 5 (t)|||2(7)J decreases monotonically in time. In the following we verify this 
fact numerically. Again, we consider the case with the initial function uo(x) = 
sin(7ra;) and the noise intensity parameter e = 1. It is not hard to prove that J(t) 
converges to zero as t — > oo. Figure [9] plots the computed J(t) as a function of t. 
The numerical result suggests that J(t) does not change anymore for t > 0.1. 

4.4. Thresholding for colored noise. In this subsection we present a compu- 
tational study of the interplay of noise and geometric evolution in ( |1.5| , which is 
beyond our theoretical results in section |3.1| and |3.2| For this purpose, we use 
driving colored noise represented by the Q- Wiener process (J € N) 



(4.2) 



W, 



j 
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FIGURE 7. Surface plots of computed solution at a fixed stochastic 
sample on the space time domains (0, 1) x (0, 0.1) (left) and (0, 1) x 
(0, 2 8 x 1CT 5 ) (right). u Q is given in p~T]) and e = V2. 




Figure 8. Surface plots of computed solution at a fixed stochastic 
sample on the space time domains (0, 1) x (0, 0.1) (left) and (0, 1) x 



(0, 2 s x 10 _J ) (right). u Q is given in (4.1 1 and e = 5 



Time from Oto 0.1 




FIGURE 9. Decay of the energy J(t) on the interval (0,0.1). 

where t > 0}j>i denotes a family of real- valued independent Wiener pro- 

cesses on (fi,.F, F, P), and ej)}^ =1 is an eigen-system of the symmetric, non- 
negative trace-class operator Q : L 2 (I) — > L 2 (I), with ej = \/2siii(jiTx). In partic- 
ular, we like to numerically address the following questions: 
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? 



(A) Thresholding: By Theorem 2.1 strong solutions of (1.5) exist for e < v2, 
and a similar result can be shown for the PDE problem with the noise (4.2 ). 
What are admissible intensities of the noise suggested by computations 
Moreover, what do the computations suggest about the stochastic MCF in 
the case of spatially white noise (i.e., qj = 1, J = oo) where no theoretical 
result is available so far? 

(B) General initial profiles: The deterministic evolution of Lipschitz initial 
graphs is well-understood. For example, the (upper) graph of two touching 
spheres may trigger non-uniqueness. What are the regularization and the 
noise excitation effects in the case of the initial data with infinite energy 
and using different noises? 

Recall that the estimate in Proposition 3.3 for V^-valued solution u S ^ n suggests 



that e > ought be sufficiently small to ensure the existence. In our test, we employ 
i 

the colored noise (4.2) with q 2 = j~ ' 6 , J = 20, and the following non-Lipschitz 



initial data: 

(4.3) u (x) = |0.5 - x\ K VxG(0,l) 



10 



where k = 0.1. In addition, we set (r,h) = (0.01,0.02) and T = |. Figure 
shows the single trajectory of the stochastic MCF plotted as graphs over the space- 
time domain with, respectively, e = 0.1,0.5, y/2. The results indicate thresholding, 
namely, the trajectories grow rapidly in time for sufficiently large values e, and the 
noise effect dominates the geometric evolution. The excitation effect of the noise on 






Figure 10. Thresholding for colored noise: Trajectories for e 
0. 1 (top left), e = 0.5 (top right), e = \/2 (bottom). 



the geometric evolution is illustrated by corresponding plots for the evolution of the 
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functional n i— > \\d x u S ^ n (oj)\\ 2 l2 vs its expectation n *— > E[\\d x u S ^ n \\ 2 L2 ] in Figure 
and [12] We observe that the geometric evolution dominates for small values o 
while the noise evolution takes over for large values of e. 





001 002 003 DM 016 016 007 003 03 1 ! 002 0D3 002 005 000 007 10 03 31 




Figure 1 1 . Geometric evolution vs colored noise evolution (g ? = 
j -0,6 , J = 20): 1st row: single trajectory for n ^ \\d x u^ n (oj)]]^ 
and e = 0.1 (left), e = 0.5 (right); 2nd row: n i-» E[||5 x ti^"||| 2 ] 
for e = 0.1 (left), e = 0.5 (right). 



4.5. Thresholding for white noise. We now consider the case of white noise 
in (3.20)-(3.21), that is, qj = 1 in < 4.2 ) and J — oo, for which the solvability of 
(1.3) is not known. Figure 13 shows the single trajectory of the stochastic MCF 
(with the same data as in section 4.4 ) plotted as graphs over the space-time domain 
with, respectively, e = 0.1, 0.5, We observe a very rapid growth of trajectories 
(numerical values range between 10 14 and 10 21 ) even for small values of e > 0. 
These numerical results suggest either a rapid growth or a finite time explosion for 
the stochastic MCF in the case of white noise. 
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